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Abstract 

We calculate the production of dileptons and photons in the presence of a nontrivial 
Polyakov loop in QCD. This is applicable to the semi-Quark Gluon Plasma (QGP), at 
temperatures above but near the critical temperature for deconfinement. The Polyakov 
loop is small in the semi-QGP, and near unity in the perturbative QGP. Working to leading 
order in the coupling constant of QGD, we find that there is a mild enhancement, ~ 20%, 
for dilepton production in the semi-QGP over that in the perturbative QGP. In contrast, we 
find that photon production is strongly suppressed in the semi-QGP, by about an order of 
magnitude, relative to the perturbative QGP. In the perturbative QGP photon production 
contains contributions from 2 —)> 2 scattering and collinear emission with the Landau- 
Pomeranchuk-Migdal (LPM) effect. In the semi-QGP we show that the two contributions 
are modified differently. The rate for 2 —?• 2 scattering is suppressed by a factor which 
depends upon the Polyakov loop. In contrast, in an SU{N) gauge theory the collinear rate 
is suppressed by 1/A^, so that the LPM effect vanishes at = oo. To leading order in the 
semi-QGP at large N, we compute the rate from 2 —)■ 2 scattering to the leading logarithmic 
order and the collinear rate to leading order. 
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I. INTRODUCTION 


In many ways, the collisions of heavy ions at high energies appear to be well 
described by thermal properties of a Quark-Gluon Plasma (QGP). Gertainly the 
bulk properties of hadrons are accurately modeled by a nearly ideal plasma, using 
hydrodynamics [1-3]. 

It is also important to consider electromagnetic probes of a QGP such as dilep¬ 
ton [4-9] and photon [10-12] production. Theoretically, these can be computed in the 
(resummed) perturbative QGP at high temperature [13], by using hadronic models 
at low temperature [14], and using the AdS/GFT correspondence [15-17]. Neither 
applies directly to Quantum GhromoDynamics (QGD) at temperatures near the the 
phase transition, at a temperature T^. 

The experiments demonstrated several phenomena which are difficult to explain 
using these methods. For dileptons, there is an enhancement at invariant masses 
below that for the p-meson. This is observed from energies at the Super Proton Syn- 
chotron (SPS) at GERN, to the Relativistic Heavy Ion Gollider (RHIG) at Brookhaven 
National Laboratory, and onto the Large Hadron Gollider (LHG) at GERN [8]. 

Another puzzle appears in the photon spectrum: there is an unexpectedly large 
elliptic flow for photons at moderate momenta, comparable to the elliptical flow 
observed for hadrons [18, 19]. This large elliptic flow for photons is very difficult to 
understand from either a perturbative analysis or from AdS/GFT. 

In this paper we consider electromagnetic signals in a matrix model of the semi- 
QGP, which is constructed to describe QGD at temperatures near and above Tc [20- 
26]. The relevant parameter is the expectation value of the Polyakov loop: properly 
normalized, the expectation value of the loop is near unity in the perturbative QGP 
[27]. Numerical simulations on the lattice [28] find that for QGD, there is no true 
phase transition, only a rapid increase in the number of degrees of freedom. For our 
purposes, whether or not there is a true phase transition is irrelevant: all that matters 
is that the (renormalized) Polyakov loop, which from the lattice is (£) ~ 0.1 at Tc, is 
small [28]. 

A brief summary of the results of this analysis has appeared previously [29] . In 
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this paper we describe the computations in full. These are straightforward, simply a 
matter of computing in the presence of a background held for the time-like component 
of the gluon vector potential, Aq. We then compute to leading order in the QCD 
coupling g. These formalisms will be explained in Sec. II. For photons, we only 
compute to leading logarithmic order, which means that we regard the logarithm of 
some large number as much larger than unity. 

In the semi-QGP, the production of colored particles is suppressed by powers 
of the Polyakov loop as T —)■ T^. This is natural, as in the pure gauge theory, there 
are no colored particles in the conhned phase. Thus one might expect that dilepton 
production is suppressed in the semi-QGP, relative to that in the perturbative phase. 
We make this comparison at the same temperature, and the same value of the QGD 
coupling, so that the ratio is only a function of the value of the Polyakov loop in 
Sec. III. In contrast to the naive expectation above, we hnd a mild enhancement of 
dilepton production in the semi-QGP, even into the conhned phase. This is because 
for an off-shell photon, it can proceed directly through a color singlet channel, a quark 
anti-quark pair. While a single quark or anti-quark is suppressed by a power of the 
Polyakov loop, a quark anti-quark pair is not. We also show that to leading order, 
a Polyakov Nambu-Jona-Lasino model [30-32] gives essentially the same result for 
dilepton production [33] as our matrix model. As we discuss, this equality is not true 
beyond leading order. 

The production of real photons, which will be analyzed in Sec. IV A, is very 
different. Kinematically, a photon on its mass shell cannot decay directly into a quark 
anti-quark pair. Therefore, the leading contribution is from a 2 to 2 scattering, which 
includes the Gompton scattering of a quark with a gluon and the pair annihilation 
of a quark anti-quark pair. These particles also can form a color singlet like the 
case of the dilepton production, but for an SU{N) gauge theory, the ratio of the 
color singlet state to the number of all the states is suppressed by 1/V^ at large N. 
Gonsequently, we hnd a strong suppression of real photon production in the semi- 
QGP. The contribution from the collinear emission of the photon, which also can 
contribute at the leading order to the photon production, is discussed in Sec. V. 
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II. SEMI-QUARK GLUON PLASMA 


A. Double line notation 


It is useful to compute the color factors using the double line basis [22], In this 
basis, as usual fundamental quarks carry a single index in the fundamental represen¬ 
tation, a = 1, • ■ ■ , iV. Gluons, however, carry a pair of fundamental indices, {ab). For 
an SU{N) group there are such pairs, and so this basis is overcomplete by one 
generator. This is compensated by introducing the operator 

K = sX - ■ ( 1 ) 

This is a projection operator. 


_ 'pafe 

' ef ’ cd • 


cd 


( 2 ) 


In the double line basis, the vertex between a quark anti-quark pair and a gluon is 
proportional to this projection operator, 

■ ( 3 ) 


The other vertices are not relevant for the present discussion. 


B. The Polyakov loop in Euclidean spacetime 

To introduce the effect of nontrivial Polyakov loop in perturbative calculation, 
we work in an effective model introduced in Ref. [22] . The Lagrangian of that model 
is the same as that in QCD with N colors, but in a mean held type approximation, 
we take the temporal component of the gluon held to be a constant, diagonal matrix, 

Af = i 5“'’ , (4) 

where g is the coupling constant. There is no background held for the spatial com¬ 
ponents of the gluon, Aj. As the gauge group is SU{N), Aq is traceless, and the sum 
of the Q’s vanishes, Q°‘ = 0. 
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The Wilson line in the temporal direction is 

f 

L(T) = V exp \ i g 

where V denotes path ordering and r is the imaginary time, r : 0 —)■ 1/T. 

To leading order in the coupling constant, a mean held approximation implies 
that we can neglect huctuations in Aq. The variable Q is naturally proportional to 
the temperature, so it is useful to introduce a dimensionless variable q, where 

= 27rT . (6) 


dr Aq{t, x) 


( 5 ) 


In this paper we shall use both the and the g“’s. For intermediate expressions 
the (5“’s are more convenient, but hnal expressions are simpler in terms of the g“’s. 
Traces of powers of the the Wilson line are Polyakov loops, 

N 

CW)^^{trL") = - (7) 

a=l 

and are gauge invariant. Since it arises frequently we write a loop without the sub¬ 
script as the hrst Polyakov loop, £ = £i. 

In general there are iV — 1 independent Q°‘A. For the problems of interest in this 
paper, though, we can perform a global color rotation to enforce that the expectation 
value of the loop ^ is real. This implies that the eigenvalues pair up as 

Q- = (-Q^, ... - 0 , .. Q^-\ gg , (s) 

where we assume that N is odd, and j = {N — l)/2. When N is even the zero 
eigenvalue is dropped, and there are j = N/2 pairs. Thus in general there are 
j independent eigenvalues. For an arbitrary value of the loop, there is no simple 
relation between these eigenvalues. 

Nevertheless, there are two exceptions. One is the perturbative QGP, where all 
g“ vanish. The other is the conhned phase of a pure gauge theory, 

AT -I- 1 —2k 

Qconf = vrT - , A; = l,...,iV. (9) 
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That is, in the conhned phase the eigenvalues are evenly distributed on the unit circle. 
The loops in the conhned phase are 




conf J 


^ n = jN- 
0 ^ nj^jN , 


( 10 ) 


for general N. This behavior is easy to understand. Loops which carry Z{N) charge 
vanish in the conhned phase of the pure gauge theory, while those which are Z{N) 
neutral do not. 

For three colors. 


= {-Q, 0, Q) = 2nT (-g, 0, q) . (11) 

The hrst Polyakov loop is then 

^ ^ (1 + 2 cos (27rg)) . (12) 

In the conhned phase of the pure gauge theory gconf = 1/3. Similarly, 

/„ = i(l + 2cos(27rng)) . (13) 

In the presence of dynamical quarks there is no rigorous dehnition of a conhned 
phase. Dynamical quarks act as a background Z{N) held, so that any Polyakov 
loop is nonzero at nonzero temperature. Nevertheless, numerical simulations on the 
lattice hnd that i is small, {£) ~ 0.1, at the phase transition, at least for three colors 
and three light havors. Thus we shall hnd it very convenient to compare results in 
the perturbative QGP to those in the conhned phase of the pure gauge theory, as a 
limiting case of how large the ehects can possibly be. 


C. Analytic continuation to Minkowski space-time 

Expanding around the background held in Eq. (4), in Euclidean spacetime the 


energy of a quark becomes 


Po ^ Po + Q°' , 

(14) 

while that of a gluon becomes 



(16) 
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where a and b are color indices of the quark and the gluon in the double line basis [22] . 
Because of the usual boundary conditions in imaginary time, the energy po for a 
fermion is an odd multiple of ttT, while that for a boson is an even multiple of vrT. 

Although the momenta for fermions and bosons are rather different in Euclidean 
spacetime, it was argued previously that the proper procedure for analytic continua¬ 
tion to Minkowski spacetime is to continue the entire Euclidean energy to —iE, where 
E is a continuous energy variable [22]. 

This has a simple but profound implication. In kinetic theory a given process is 
given by an integral over phase space of the square of a matrix elements times products 
of statistical distribution functions. Since the energies in Minkowski spacetime are as 
usual, then, for processes in which all the momenta are hard, the only change is in the 
Q-dependence of the statistical distribution functions. For processes involving soft 
momenta, it is also necessary to include the Q-dependence of the hard thermal loops as 
well [22] . We shall illustrate these general expectations by our explicit computations. 
It also suggests that it may be useful to treat the semi-QGP in kinetic theory, as for 
the perturbative QGP [34, 35]. 

The background gluon held acts as an imaginary chemical potential for colored 
particle, so that the statistical distribution functions for the quark, anti-quark, and 
gluon are, respectively. 


na{E) = 
nab{E) = 


^{E-iQ‘^)/T I 

1 

g(E-i(Q“-Qh)/'r 


, na{E) = 


g{E+iQ-)/T X ’ 


(16) 


Notice that the sign of Q for the anti-quark, +iQ°‘ in na{E), is opposite to that for 
the quark, —in na{E). This is just like the change in sign for a quark chemical 
potential which is real. When the Q°' = 0, of course these reduce to the usual Fermi- 
Dirac and Bose-Einstein distribution functions. 

For future reference, it is useful to compute the statistical distribution functions, 
summed over all colors, in the conhned phase of a pure gauge theory, Eqs. (9) and 
(10). For the quark distribution function, this is 

N 




TV ^ + 1 + 1 ’ 

a=l 


( 17 ) 







while that of the gluon distribution function is 



a,6=1 


1 

e(®“**3conf+*Qconf)/'^ — 


1 

^NE/T _ I • 


( 18 ) 


In the conhned phase of the pure gauge theory, the only loops which contribute are 
those which wrap around a multiple of N times. These can be considered as a type 
of “baryon”, albeit in the pure gauge theory. Consequently, the energy which enters 
in the right hand side of Eqs. (17) and (18) is not E, but N times E. This rescaling 
of the energy will be seen to help explain the suppression of photon production at 
large iV, Eq. (125). 


D. Relation to lattice results 


Strictly speaking, and thus should be determined dynamically from our 
model itself. Instead, in this paper, we determine these quantities from numerical 
simulations on the lattice, following Ref. [26]. First, in order to extract the Polyakov 
loop determined by nonperturbative dynamics, it is necessary to remove perturbative 
corrections from the expectation value of the loop [27], 


£(Q = 0) = 1 + 6i{Q = 0) , 


(19) 


-^n2 + iV n — + - 

2 I T 4 


+ 0(9= 


where Cf = (iV^ — 1)/(2A^) is the Casimir for the fundamental representation, mo is 
the Debye mass of the gluon, and Nf is number of quark flavors. We use the running 
coupling constant calculated in the modihed minimal subtraction scheme at two-loop 
order, and the expression of the Debye mass at one-loop order [27]: 

n "1 


C/2 = 247r2 


{IIN - 2Nf) jin 


47rT 

^MS 


-1e\ + Nf{A In 2 - 1) - 


lliV 


= 


{2N + Nf)ATi'^T‘^ (lliV - 2Nf) | In ( ^ 

L I V ^MS / 

5iV2 + iv| + 9Nf/{2N) 


+ 4iV/ln2 


2N + Nf 


( 20 ) 


Here is the renormalization mass scale, in the modihed minimal subtraction 
scheme, and je — 0.57721 is Euler’s constant. 
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FIG. 1. Left panel; The Polyakov loop {1) determined from the lattice calculation [28], the 
Polyakov loop in which the perturbative correction is removed (Iq) as a function of T. Right 
panel: Q as a function of T. We set = Tc/l-Sh, where Tc = 170 MeV. 


Equation (19) shows that a hnite renormalization gives £{Q = 0) > 1. We 
assume that perturbative corrections exponentiate, 

i(Q) = • (21) 

We take i from numerical simulations of lattice QCD [28], and calculate io from 
Eq. (21), to obtain Q from Eq. (12). These quantities are plotted in Fig. (1), by 
setting Tlj^ = Tc/1.35. We see that io is different from unity even around ~ 3Tc, 
where Tc ~ 170 MeV is the pseudo-critical temperature of the phase transition [28]. 


III. DILEPTON PRODUCTION RATE 


A. Computation to leading order 


We calculate the production rate of dileptons when 7 ^ 0 in this subsection 
To leading order in Oem, the production rate is 

dr 


W'T(P), 


( 22 ) 


(i-ip 247r‘p2 

where LF^i,(P) is the Wightman correlator for two electromagnetic currents, 

W,,iP) = [ d^x ( jM Ux) ) , (23) 
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where = e 'Yhf'4^f i with 'ip is the quark operator with flavor index /. In thermal 
equilibrium, W^i, is related to the imaginary part of the retarded photon self-energy 
as 


W,JP) = - 2 n(E) Im n^,{P) , 


(24) 


with 


= -' j 9(^°) ( [jj.(a:).i.-(0)] > • 


(25) 


Here P = Pi -|- P 2 with Pi and P 2 being the momenta of the two leptons. 

At the leading order in the QCD coupling constant g, the contribution is ob¬ 
tained by 1 —)■ 2 processes, illustrated in Fig. (2). In this process, a quark anti-quark 
pair becomes a virtual photon, which then decays to a dilepton pair. This gives 


dr 

'¥p 




247r4p2 


E 

/.spin 

N 


d?ki 1 r d?k2 1 

(27r)3 2Pi j (27r)3 2 P 2 


(27r)"5W(P-Pi-P2) 


X \M\^ na{E,)ME2 


a=l 


(26) 


where / is a subscript for flavor running from 1 to iVj. We use the spacetime signature 

(-1- —-) in this paper; four-momenta are denoted by capital letters, P^ = {E,p), 

K'p = (Pi, ki), K 2 = (P 2 , ^ 2 )- The quark anti-quark pair is produced on it mass shell, 
Kf = = 0, so Pi = \ki\ = ki, P 2 = 1 ^ 2 ! = ^ 2 , and P is time-like, P^ > 0. Without 

loss of generality we can assume that the (virtual) photon energy is positive, P > 0. 
Here na(Pi) and na{E 2 ) are the statistical distribution functions for the quark and 
anti-quark in Eq. (16). The square of matrix element is 


^|Wl|2 = 8e2gJiFi-P2=4e2gJp2, (27) 


where we have used = 0, and O'/ is the electromagnetic charge of the quark 

with flavor / in the unit of e. 

The result when Q = 0 is well known [9]: 

= ^ ^ • (28) 
/ 


dr 

#p 
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For three flavors of quarks, YIjQ} — 2/3. Here 


h{E,p) = 1 - 


2T, 
— In 
P 


n+e-P-/^\ 

{l + e-P+/^) ’ 


and 


(29) 


P± = ^{E±p) 


(30) 


is the range of the quark momenta. Especially, when the dilepton pair is produced 
at rest, p = 0, the quark anti-quark pair are produced back to back, with ki = —k 2 - 
Their energies are equal, Ei = E 2 = E/2, and there is no integral over the quark 
momentum. The expression then reduces to 


dr 


d^P 


a: 


Q=o 


127r4 


Nn^(E/2). 

f 


(31) 


This is natural, as the product of a Fermi-Dirac distribution function for the quark 
and anti-quark appears. 

Equation (26) illustrates our comment in Sec. II C, that for hard momenta the 
only change when 7 ^ 0 is in the change in the statistical distribution functions. To 
handle the Q-dependence of the u’s it is useful to note that 


naiEi) UaiE^) = n{E) (1 - na{E,) - na{E 2 ]) , 


(32) 


remembering that E = Ei + E 2 . This identity is familiar from when Q°' = 0. 

Using this, we can derive 

/ ^2 _ 2 \ ^ rp+ 

Im = ttem q} - - - j YlJ dk{l- na{k) - na{E - k)) . (33) 

Here we wrote the photon retarded self-energy instead of the dilepton production rate 
for future convenience. We use energy-momentum conservation to write Ei = k and 
E 2 = E — k. This expression is useful when we compare to the results of Ref. [33] at 
the end of this section, see Eq. (44). 

To leading order, we can write the rate for dilepton production when 7 ^ 0 
as a momentum dependent factor times that for = 0 , 


dr 


Q70 


fi-M 


dr 


d^p 


Q=0 


(34) 
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FIG. 2. The 1 to 2 process which results the production of dilepton. The solid line denotes 
the quark while the wavy line denotes the (virtual) photon. 


where 


KE,V) fiAQ) 


1 A/ 2T /I+ e-(P— 

N ~ y I l + e-(P+-*Q“)/^ / / 

a=l \ J \ / / 


(35) 


This result can be evaluated by expanding in powers of exp(—(p:p — iQ°‘)/T), and 
performing the sum over a to obtain a series of Polyakov loops. For general N 
all independent Polyakov loops, which run from £i to enter. The resulting 

expression is not very illuminating. 

There are two cases in which one can obtain simple results. One is the conhned 
phase of the pure gauge theory, Eqs. (9) and (10). Then only loops which are a 
multiple of N contribute, so that 

h{E,p) faiQconf) = 1 - — In j • (36) 


Another special case is three colors. Then one can rewrite fij so that only the 
hrst Polyakov loop appears. 


h{E,p) fiiiQ) 


2 T /1 + 3 £ + 3i e-^P-/^ + e-3p-/^\ 

3p yl + 3ie~P+/'P + 3ie~‘^P+/'^ + 


(37) 


Of course fij{0) = 1 in the perturbative QGP, when £ = 1. In the conhned phase 
where £ = 0, this agrees with the result in Eq. (36). 

In Fig. (3) we plot fijiQ) as a function of temperature for three colors. We do 
this for back to back dileptons, p = 0, with E = 1 GeV. We see that the production 
of dileptons is not suppressed by the effect of the Polyakov loop, but moderately 
enhanced, by ~ 20%, at low temperatures T ~ 300 MeV in the semi-QGP. 
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FIG. 3. The ratio of dilepton production in the semi-QGP versus the perturbative QGP, 
ffi in Eq. (37), as a function of temperature. The dileptons are back to back, p = 0, with a 
total energy E = 1 GeV. 


This enhancement is rather unexpected. While the probability to produce either 
a single quark or anti-quark is small when the loop is small, that to produce a quark 
anti-quark pair is greater in the semi-QGP than the perturbative QGP. 


B. Enhancment of dilepton production in the confined phase versus the 
perturbative Quark-Gluon Plasma 

To better understand the enhancement of dilepton production in the semi-QGP, 
relative to that in the perturbative QGP, we consider dilepton production for infinite 
N, comparing the confined phase to the perturbative QGP. 

To simplify the analysis we consider dileptons which are produced back to back. 
This is most useful, because if the total spatial momentum of the pair vanishes, p = 0, 
then each dilepton carries the same energy, E/2, and we can ignore the integral over 
phase space as a common factor, independent of the The effects of confinement, 
represented by 7 ^ 0 , can then be included just by computing the sum over the 
statistical distribution functions in Eq. (26), 

N ^11 

^ na{E/2) na{E/2) — ^ ^{El2-iQ‘^)IT I g(E/2+iQ“)/T ]_ ' 
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We note that, since the background held acts like a chemical potential for color, albeit 
imaginary, the sign for is opposite between the quark and the anti-quark. 

In the perturbative phase, = 0, and Eq. (38) is just = iVn(E/2)^, as appears 
in Eq. (31). 

In the semi-QGP, Eq. (38) is computed by expanding each statistical distribu¬ 
tion function in powers of exp((E/2 =f iQ°‘)/T), 

N OO OO 

EEE Qxp {—{{m + m')E/2 + i{m — m')Q°‘)/T) . (39) 

a=l m=l m'=l 

This sum is especially easy to compute in the conhned phase at inhnite N. In that 
case, if m 7 ^ m' the sum over a gives ^exp(i(m — m')Q°'/T)] this is the Polyakov 
loop whose contribution vanishes at large N. The only nonzero contributions 

are from terms where m = m'. For the terms in Eq. (39) where m = m', though, 
the dependence on drops out, cancelling identically between the quark and the 
anti-quark. The sums over a and m are then independent, and easy to do, 

E E N . («) 

a=l m=l 

which does not vanish at large N. That is, while we start with only Fermi-Dirac 
distribution functions with ^ 0, in the conhned phase at inhnite N we end up with 
a Bose-Einstein distribution function, which corresponds to the mesonic distribution 
function instead of the quark and the anti-quark. We also note that, previously we 
showed that the cancellation of the phases of the quark and the anti-quark are essential 
for the non-suppression of the dilepton rate at large N by using the Boltzmann 
approximation [29]. From the discussion above, we see that the cancellation (m = m') 
is important also in the case that we do not use the Boltzmann approximation. 

This is a type of statistical conhnement. Our simple model does not have true 
bound states, but there is a remnant of a bound state from the statistical sum over 
the Q“’s. It is this sum in Eq. (39) which generates the Bose-Einstein distribution 
function in Eq. (40). 

Thus in the conhned phase at inhnite N, 

?7.(E/) 

flliQconf)N=oo = 
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We note that this result also can be obtained by taking p —)■ 0 limit in Eq. (35). This 
demonstrates a few interesting features. First, fijiQconi) n=oo is always larger than 
unity. Second, at low energy, the Bose-Einstein distribution function is enhanced, 
n{E) ~ T/E, while the Fermi-Dirac distribution function is constant, so 

4T 

ffliQ coni)N=oo = , E <^T . (42) 

Thus under the given assumptions, at small energies dilepton production in the con- 
hned phase dominates that from the perturbative Quark-Gluon Plasma. This occurs 
because statistical conhnement generates conhned “bosons” from quark anti-quark 
pairs, and these conhned bosons become over-occupied when their energies are much 
smaller than temperature. This occurs even though the probability to produce a single 
quark, or anti-quark, is strictly zero in the conhned phase at inhnite N. Nevertheless, 
we note that, when E < gT, we need to calculate with the HTL resummation [36] 
instead of our calculation, so our result Eq. (42) can be altered in that energy region. 

More generally, that the ratio of dilepton production in the conhned phase to 
that in the perturbative QGP, fijiQ), is of order one, indicates that at all temperature 
dilepton production is of order N. This is one example of quark-hadron duality [37]. 

A similar enhancement of dilepton production in the conhned phase was found 
previously by Lee, Wirstram, Zahed, and Hansson [6]. They considered a condensate 
for Aq. We can take our result in Eq. (35), and expand up to quadratic order in the 
Q“’s, to obtain 

1 ^ 

HE,p) fiiiQ) Si h{E,p) + («(P-)(1 -n(p-)) -«(p+)(1 -n(p+))) . 

( 43 ) 

in agreement with Eq. (7) of Ref. [6]. These authors suggested that the enhancement 
of dilepton production in the conhned phase may be related to the excess seen in 
heavy ion collisions for dilepton masses below that for the p-meson [8]. 

We can also make contact with results from Polyakov Nambu-Jona-Lasino 
(PNJL) models [30-32], especially with the computation of dilepton production by 
Islam, Majumder, Haque, and Mustafa [33]. To do so we need a simple identity. For 
three colors the sum of the Fermi-Dirac distribution function, with the and i as 
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in Eqs. (11) and (12), obeys 


( 44 ) 



1 + 3 £ + 3 £ e-2®/'^ + ' 


i e + 2 £ e + e 


a=l 


a=l 


In the PNJL models of Refs. [30-33], when I ^ \ the effective statistical dis¬ 
tribution function is defined as the right hand side of Eq. (44); e.g., Eqs. (67) and 
( 68 ) of Ref. [31]. In Refs. [31-33] this effective distribution function was obtained by 
taking the derivative of the free energy, when i 1, with respect to a given energy 
E. Since the free energy involves a sum over all colors, it is clear that dehning the 
effective statistical distribution function in this way automatically gives a sum over 
all na{E) (or na{E)), which appears on the left hand side of Eq. (44). (The identity 
of Eq. (44) holds for the case of zero quark chemical potential. Then we can dehne 
i to be real, and Yha^o-^E) and equal. At nonzero quark chemical 

potential the loops in the fundamental and anti-fundamental representations are not 
equal, £3 7 ^ £3 [38]. In this instance, identities similar to Eq. (44) hold for J2a^a{E) 
and J2a^a{.E) separately, and are again equal to those in the PNJL model, [31-33].) 

In our matrix model the sum over the statistical distribution functions with all 
( 5 “, J2a^a,{E) and '^a^a{E), enters naturally when we sum over all quark colors, 
Eqs. (32) and (35). In the PNJL model calculation done in Ref. [33], Imi7(^ is given 
in Eq. (4.36), and by Eq. (4.46). Taking the quarks to be massless, and using 

the fact that dp na{E — p) = dp na{E), it can be shown that their result for 
Imll^^ coincides identically with our Eq. (33). 

We emphasize that the equality between our results and the PNJL model [33] 
is valid only to leading order. In both cases, at leading order dilepton production is 
only a function of the Polyakov loop and the temperature. (As well as the quark mass 
and chemical potential, if one chooses to add them.) The results will certainly differ 
beyond leading order, and depend strongly upon the details of each effective model. 

We note that, results for dilepton production at nonzero quark masses were 
computed in Ref. [33] and by Satow and Weise [39] . There is a relatively mild depen¬ 
dence on the quark masses, apart from obvious kinematical constraints, such as the 
energy of the photon has to be greater than twice the quark mass. 
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IV. PHOTON PRODUCTION 


A. Overview 


To leading order in cTem, the photon rate in the QGP is 


dF-y 


1 

2(27r)3 




( 45 ) 


Since a photon on its mass shell cannot decay directly to a quark anti-quark pair, 
this quantity vanishes at one-loop order. 

In our model the hrst nonzero contribution occurs at two-loop order, from the 
diagrams shown in Fig. (4). Cutting the diagrams we obtain 2 —)■ 2 processes, which 
are Compton scattering and pair annihilation, both of order e^g^. We note that, 
consequently, and unlike the case of dilepton production to leading order, the results 
which we hnd have no direct correspondence in a PNJL model. One could compute 
photon production in a PNJL model, but since these models do not have dynamical 
gluons, the results will be very different from our matrix model. Each of these two 
processes has an infrared divergence when the momentum exchanged becomes soft [4] . 
The divergence is removed by using a resummed quark propagator for soft momenta, 
corresponding to the uncutted lower quark line in the left diagram of Fig. (4), for 
example. 

It was later realized that there exists another kinematic regime which con¬ 
tributes at the same order [5]. This corresponds to the case when the photon becomes 
collinear with quarks in the loop in the two diagrams of Fig. (4), more precisely, when 
the longitudinal momenta (dehned with respect to photon momentum) of quarks re¬ 
main hard, ~ T, and the transverse momenta are soft, ~ gT. Despite the reduced 
phase space, due to collinear enhancement, this regime was found to contribute equally 
as the 2 —)■ 2 processes, in the analysis by Arnold, Moore and Yaffe (AMY) [5]: The 
collinear regime in this diagram also gives an overall e^g^ contribution to the photon 
emission rate. To clarify terminology, we will refer 2 -f-)- 2 rate as the contribution 
from Fig. (4), excluding the rate in the collinear regime. We refer to the rest of the 
contributions ~ e^g"^ as the collinear rate. 

It turns out that the collinear rate goes beyond two-loop order: i.e. further 
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FIG. 4. Two loops diagrams contributing to photon self-energy. 

additional soft gluons ladders in the collinear regime still contributes at the same or¬ 
der. Thus, the story is further complicated by interference among different diagrams. 
Physically, it is because that, the formation time of a photon, tp ~ l/{g‘^T), is com¬ 
parable to the mean free path for quarks, A ~ l/{g‘^T). Since these two scales are 
similar, interference effects between scattering with multiple gluons must be included, 
which is the Landau-Pomeranchuk-Migdal (LPM) effect [5]. Different diagrams add 
destructively, so that the LPM effect leads to additional suppression of collinear pho¬ 
ton rate by at large photon momentum p ^ T. 

In this section we compute the production of real photons with large momentum 
in the presence of a nontrivial Polyakov loop. We begin by reviewing the computation 
of photon production to leading order in perturbation theory for 2 —)■ 2 processes. 
We then generalize this to ^ 0. In contrast to dilepton production, we hnd that 
photon production is strongly suppressed in the conhned phase. We give a simple 
explanation for this in terms of the initial state of the scattering. 

We then give a detailed computation of the leading contributions to the collinear 
rate when Q°‘ ^ 0. In the presence of a nontrivial loop, the thermal mass of the quark 
is suppressed by a loop dependent factor, but it remains ~ gy/NT (here we explicitly 
wrote the N dependence in the large N limit). In contrast, the damping rate is 
suppressed by a factor of 1/N. Consequently, the mean free path of a quark or gluon 
is much larger, A ~ l/{g'^T), not l/{g'^NT) as in the Qa = 0 case. This implies that 
the LPM effect can be neglected at large N. 

We compute the collinear processes when 7 ^ 0 at large N. Doing so, we hnd 
that for three colors, the result is not that small, at least for physically reasonable 
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values of the QCD coupling constant. Nevertheless, we hnd the result illuminating, 
to show how results can change in the semi-QGP. 


B. Hard momentum exchange with trivial Polyakov loop 


To establish notation on kinematics, we hrst review the computation of the 
differential photon rate for 2 -f-)- 2 processes at hard momentum exchange, in the case 
of Q°' = 0 [4]. In kinetic theory, this is given by 


dr 



d^ki d^k2d^P' 
{2nf9,EiE2 




P-P') 5{{P'f)d{E') 


X n{Ei) n{E 2 ) {l±n{E')) \M\l . 


(46) 


The summation i represents the contribution of Compton scattering and pair annihi¬ 
lation, whose diagrams are shown in Fig. (5). The statistical factors n{Ei), ''^(-^ 2 ), and 
n{E') can refer to either Fermi-Dirac or Bose-Einstein factors, depending upon the 
particular process. For Compton scattering, the statistical factor above is 1 — n{E'), 
which corresponds to Pauli blocking; for pair annihilation, the corresponding factor 
1 -|- n{E'), which represents Bose enhancement. 

The incoming momenta are Ki = {Ei,ki) and K 2 = (^ 2 ,^ 2 ) the outgoing 
mometum P' = {E',p'), and P = {E,p) is the photon momentum. We assume all 
particles are massless, so Ei = \ki\, etc. Whether the incoming or outgoing momenta 
are quarks or gluons depends upon the process considered. In this paper, we consider 
the case that the photon energy is much larger than temperature, E 3 > T. 

It is convenient to introduce the Mandelstam variables, 

S={K, + K2)\ 

t = (iPi - P)2 , (47) 

U={K2- Pf . 

With our kinematics, 

s > 0 ; t ,u < 0 . (48) 
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We decompose the incoming momenta ki and ^2 into components parallel and per- 
pendicnlar to the photon momentnm p, with 


7 II t /T I \2 t^l 

= 7^ + El ^ (^1 ) = -> 

2p 4p^ p 

7II u „ /, I x9 uE^ 

k\ = - + E^, [kif = - - 


2p 


4p2 


p 


(49) 


and where ki ■ ^2 = k^k^ cos{(j)i — 02)- We can then convert the variables of 

integration as 

1 \ E E 

(fki d^k 2 = - d(j)i d02 dkl dkl dik^Y = - d0i d02 \ ^ dt du dEi dE 2 . (50) 

4 4 

The integrand only depends on 0 = 0i — 02 throngh 6{{P'Y)- 


6{P'j = 6{2 EiE 2 — 2k\k'2 — cos0 — s) . 


(51) 


The angular integrals are easily done to give the following result: 


r d^kid^k 2 

J SEIE 2 



dtdudxdy {277) 1 

32p -\-hy + c ’ 


where we dehne 


X = El + E2 ; y = El- E2 , 


and 



b={^-p){t^-uY, 

^ /j. \2 2 I 2 2 2 j. 

c = —-{t — u) X + ps X — p S — uts . 


(52) 


(53) 


(54) 


We start with the integral over y. Let y± be the solutions of the quadratic form 
in y, ay\ + hy±Pc = 0. The integral over y runs from y_ to p+, where ay'^ + hy + c > 0. 

In considering the quadratic form in y, we assumed that 6^ — 4ac > 0. A bit 
of algebra shows that this determines the range for x to be x > p + s/(4p). As the 
energy E' = Ei + E 2 — p = x — p, we automatically satisfy the condition that this 
particle has positive energy, E' > 0, and can set 0{E') = 1 in Eq. (46). 
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Since we assume that the incoming momenta are hard, the distribution func¬ 
tions, n{Ei) and n{E 2 ), can be replaced by their Boltzmann forms, exp{—Ei/T) 
and exp{—E 2 /T). Consequently, the product of statistical distribution functions in 
Eq. (46) reduces to 

n{EME,)(l ± n{E')) ~ (l ± ^ 577 ^) = (l ± 777 = 75 ) 5 ^ 

(55) 



This vastly simplihes the integral over phase space. In general, the product in Eq. (55) 
is a function of both sum and difference of the energies, x and y. For hard momenta, 
though, this reduces just to a function of the sum, of x. In appendix A, we show that 
corrections to Eq. (55) are in fact exponentially suppressed, as one would expect. 

This allows us to immediately perform the integral over y. Although the coef- 
hcients b and c, and y±, are all functions of x, in the end we obtain simply 


ry+ 


dy 


y- y/ay"^ + by+ c 


sm 


2ay + b 
— 4ac 


y+ 


TT 


(56) 


We can then readily evaluate the integral over x. 


dx e (^1 ± ^ T In (1 T e . (57) 

Therefore, the phase space integrals give 

= j dtdu In (1 T . (58) 

To proceed, we consider Compton scattering and pair annihilation separately, 
since it involves a calculation of the matrix element squared. For Compton scattering 
off of quarks and antiquarks, the squared amplitude is given by 

|A1p = 2 (47r)V„ a, (-8) . (59) 


The hrst term in Eq. (59), ~ s/t, is logarithmically divergent when integrated over 
t. The second term, ^ t/s, does not produce a logarithmic divergence. As we show 
below, it is the logarithmic divergence that gives rise to leading logarithmic results 
in photon production, and we can ignore the second term. 
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Remembering that t is negative, the logarithmic divergence happens for small 
—t, and invalidates the kinetic theory description. The standard treatment is to 
introduce an IR cutoff /i for the spatial component of the exchanged momentum, 
1^1 ~ P I > h- We assume that this cutoff lies between the hard and soft scales in the 
problem, p S> gT and p ^ T. Near zero, the integral over t is modihed as follows: 

1^1-P I >P ^ t < (-El - . (60) 


On the other hand, the integration range of y is given by + &?/ + c > 0, which as 
f —)■ 0 takes the following form 

2\Ei — p\ = \x + y — 2p\ <2 \j-{Ap^ + s — Apx) . (61) 

Comparing Eqs. (60) and (61) and noting x > p, the lower cutoff on —t is 


-t > 


s 

4p(a; — p) 



(62) 


Since we compute only to leading logarithmic accuracy, in the integral over —t we 
can simply take the lower limit to be /r^, to obtain 

This leaves an integral over u. However, since s = —t — u, we can trade this for an 
integral over s. The final s-integral becomes 

-^ = -^ (4pT) In > (64) 

where we have replaced ln(s//i^) by ln(pT//i^) and extend the lower bound of the 
integration to zero. This is justihed as to leading logarithmic order the region of 
integration is s ~ pT S> 

Collecting everything together, we obtain 


dr 




Q^em Cl* / , x2 
487r2 ^ 


1) In 
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The case of annihilation proceeds similarly. The sqnared amplitnde is given by 


]\j2 _ 1 

l-MP = W^aemOis - - - 

/ 



( 66 ) 


Since the integrand is symmetric in t and u, both t and n-channels contribnte the 
same to leading logarithmic order. The integral in the t-channel becomes 


u 


d{-t) - = h 


s + s In — . 


(67) 


We again keep only the logarithm and nse the same trick as in Eq. (64) to obtain 
the leading logarithmic resnlt. Note that there is Bose-Einstein enhancement for the 
annihilation process: 


r y 1" (1 - i°( A) 

J ^2 S H 

^ (^) L ~ = -y (4pT) In . (68) 

The M-channel gives an identical contribntion. 

Collecting everything together, the combination of Compton scattering in the 
t channel, and pair annhilation in the t and u channels, is 


dr 




Q^em Ois / j,j-2 

IbTT^ 


1 ) 



(69) 


C. Hard momentum exchange with nontrivial Polyakov loop 

In the previons section, we compnted the matrix elements for the diagrams 
which contribnte to photon prodnction at leading logarithmic order. Once we work 
in terms of Minkowski variables, there is no change in compnting in the presence of 
a backgronnd held for the Polyakov loop. 

The only change in a backgronnd held arises from the modihcation of the sta¬ 
tistical distribntion fnnctions. We start with the case of Compton scattering, as 
illustrated in the hgure in the left hand side of Fig. (6). In this case, the incoming 
momenta are those of a gluon, with momentum Ki, and a quark, with momentum K 2 . 
Consequently, in the statistical distribution functions we replace the gluon energy as 
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FIG. 5. The diagrams for the Compton scattering (left) and the pair annihilation (right). 
The solid line corresponds to a quark, the wavy line to the photon, and the curly line to 
the gluon, respectively. The Compton scattering includes s and t channel processes while 
pair annihilation includes the t and u channel processes. 




FIG. 6. The color labeling of Compton scattering (left) and pair annihilation (right). The 
double line corresponds to a gluon, a single line to a quark or antiquark. The color flow does 
not involve color neutral photon, which we still indicate with a dashed line. The quark- 
gluon vertices are drawn as an empty circle. They have their own graphic representation 
[22], but this is not needed here. The Feynman diagrams obtained by crossing symmetry 
are identical in the flow of colors. 

El ^ El + i{Q^ — Q^), while the quark energy E 2 ^ E 2 + iQ^. Similarly, the energy 
of the outgoing quark becomes E' ^ E' + iQ^. 

With the color labeling in Fig. (6), the thermal distribution functions when 
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7^ 0 are 


dx I X_ 


>p+s/{Ap) 


Q{x-p+iQe)IT _|_ ^ 


— g—*(Qi)+Qc—Qd)/T' 


poo ^ 

/ dx I 
Jp+s/ip 


-x/T ^_^p^-n{x-p+iQe)/T 


_ _ { _1 '\n+l 

— g-i{Qb+Qc-Qd)/T ^ ' J' q-p/T ^-nsl{ipT)^-i{n-l)QelT 

^ n 


n=l 


To obtain the leading logarithmic result, we recall Eq. (64): Eq. (70) should be 
integrated over s. The s-dependent factor gives rise to an additional factor of 1/n: 


ds exp 


n \ 4pT 


(71) 


ApT j n 

It is sufficient to calculate the ratio of photon rate with 7^ 0 to that in the 
perturbative limit, (5“ = 0. We will thus only keep track of Q-dependent factor 
^^(-l)”+Wxp(-i(Qb + Qc - Qd+ {n - l)Qe))/T)/n^. To proceed, we then need 
the form of the quark-gluon vertex in the double line notation [22], appearing in the 
matrix element squared, Eq. (3). We then multiply Eq. (70) by the product of two 
quark-gluon vertices. 


b,c,d,e 


- 1 ) 


n+l 


^ — i{Qb-\-Qc~Qd~^{^~^)Qe)/'^ 


n=l 




1 


^ ^ I ^bd ^ce ^bd ^ce ^cd ^be “ 1 “ j^2 ^ ^ 


b^c^d^e 

Ar2_i 
2N 


E 

n=l 


N 

(_l)n+l 


(_l)n+l 


n=l 




tr 17' 






(72) 


When all Q's are zero, this reduces to 


(- 1 ) 


n+l 


n=l 




N‘^-l\ + 
12 


(73) 


For Compton scattering, the ratio of this contribution when 7^ 0, to that for 
= 0, is just the ratio of Eqs. (72) and (73), 


12 


/comp(Q) - — 

vr^ ^—' 


(- 1 ) 


n+l 


n=l 




( 74 ) 


where + is the n-th Polyakov loop in Eq. (7). 
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The case of annihilation is similar. In the presence of a backgronnd color charge, 
the thermal distribution becomes 


dx ( 1 + 


lp+s/{4p) 


Q{x-p+iQc-iQd)IT _ 


-p/T „-ns/{4:pT) -i{Qi,-Qs+{n-l){Qc-Qd))/T 


= V- Te e -e 

^ n 


n=l 


(75) 


Again, the integration of e over s picks up an additional factor of 1/n. The 

color sum for scattering in the f-channel becomes 


b,c,d,e n=l 



iiQb-Qe + irt-l)iQc-Qd))/T 


1 

2 


^ ^ ^ ^bc ^de jY ^bc ^de ^cd ^be T 
b^c^d^e ^ 



^-i(Qb-Q<L+{n.-^){Qc-Qd))/T 


1 

2 


OO -| 

n=l 


(76) 


When all Q's are zero, Eq. (76) becomes 



(77) 


Scattering in the u-channel gives a result identical to that in the f-channel. Therefore, 
the suppression factor for annihilation is given by 


/pair(Q) 


iV2 




TT^ 


n=l 


(78) 


Remember that Compton scattering is 1/3 of the total for 2 —)• 2 scatterings, Eqs. (65) 
and (69). Summing over Compoton scattering and pair annihilation, to leading loga¬ 
rithmic order, we obtain the contribution from 2 —)■ 2 scattering from hard momenta 
in the semi-QGP, 

^ ^ I"" (^) A(Q) > 

fl{Q) = ^ (/comp(g) + 2/pair(g)) • 
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These expressions can be more simply expressed when = 3 in terms of (5“ 
2vrT(-g,0,g), Eq. (13): 


/comp(?) 1 8^5 

/pair(g) = l-6g + 9g^, 

f-fiQ) = 1 - 4g+ y . 


(80) 


The resnlts for more than three colors are similar, simple qnadratic polynomials in 
the Q“’s. That for /comp(Q) involves the while that for /pair(Q) is a fnnction of 
the differences, — Q^. 

We also note that exactly the same fnnctions of g“ enter into collisional energy 
loss for a heavy qnark in the semi-QGP. Becanse of the historical convention, the 
fnnction for Compton scattering in photon prodnction, /comp(g), is identical to that 
for Conlomb scattering of a heavy qnark, Eq. (33) of Ref. [26]. Similarly, the fnnc¬ 
tion for pair annhilation in photon prodnction, /pair(Q), is the same fnnction as for 
Compton scattering of a heavy qnark, Eq. (45) of Ref. [26]. While these two fnnctions 
are the same, in detail they enter differently into collisional energy loss for a heavy 
qnark, times different logarithms of the energy. 

D. Soft momentum exchange 

We now compnte the contribntion to photon prodnction when the momentnm 
exchanged is soft. This case is simpler than when the momentnm exchanged is hard, 
and so we treat the case of a nontrivial Polyakov loop at the ontset. 

We follow the analysis of Baier, Nakkagawa, Niegawa, and Redlich [4]. We 
begin the compntation in imaginary time, and then analytically continne the external 
momentnm. The photon self-energy in the imaginary time is 



(81) 


The overall factor of two arises because K or K — P can be a soft momentum: we 
have chosen only K to be soft. Thus the momenta K — P is hard, so we can use the 
bare quark propagator, Sa{K — P). For the quark with soft momenta it is necessary 
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to use a propagator, Sl{K), which is resummed with Hard Thermal Loops (HTLs) 
in the presence of 7 ^ 0 [ 22 ], 


Sa{K) 

s:{K) 



duj p{oj, k) 
2tt oj — ik^ 
duj p‘Joj,k) 
uj — ikf) 


(82) 

(83) 


with ko = ko + Qa and the quark spectral functions are 


p{uj, k) = , 

*/ A * / A 7° + 7'^ * / 

P«Afc) = -;;-P„+(w.i) + - -;;- Pa-(‘^,k) , 


(84) 

(85) 


2 ' 2 
where eiu) is the sign function. We note that the bare quark spectral function p(a;, k) 
does not have its color index. The HTL spectral functions are a sum of pole and cut 
terms, 

p*±(a;, k) = 27r [Z±a{k) 5{u: - u:±a{k)) + Z:^a{k) 5(n; + uj^a{k))] 

+ B{e - k) . 

The quark quasi-particles have a thermal mass rriqka, a dispersion relation u)±a{k), 
and residue Z±a{k) {k = |fc|). Explicitly, 

uUk) - e 


( 86 ) 


Z±a{k) = 


2ml^a 


(^±a{k) Tk = 


m 


qk< 2 - 


lT^)Qof^l±l 


where 


k 




k 


(87) 

( 88 ) 

(89) 


The explicit form of the cut term from Landau damping, is irrelevant for 

our analysis. The result for the quark quasi-particle mass is given later. 
Introducing a spectral representation for the propagators, 

d?k r dui r d(jj2 

^ 7 / r27ri3 / ^ ^ 

(nZi,)k\ —'n.h.ii'l') 

7Va(^n ^)7/4P(a;2, A; - ^ 






+ ui- U 2 


(90) 
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Since for massless quarks their spectral density has only a vector component, 


tr 


7Va(^i>^)7/.P(^2,/c -p) 


= - 2 tr 


P*a{^i:k)p{u:2,k - ^ 


(91) 


Now we compute the discontinuity in the amplitude, as we analytically continue 
the photon energy po —)■ —iE ± e, for inhnitesimal e. 




5 ^ (n„(w - E) - 7ia(w)) 

ZTT 


X tr 




(92) 


where Disci7(;‘(P) = [i7^(P + ie,p) — — ze,p)]/(2i). Since the photon is a singlet 

under color, there is no ambiguity in how we do the analytic continuation for the 
photon energy. We have also used the fact that the spectral function is real. When 
the Q“’s vanish, this discontinuity is the same as the imaginary part of the retarded 
self-energy. When the 7 ^ 0, however, if we were to compute the imaginary part, we 
would also obtain contributions from the imaginary parts of the statistical distribution 
functions, which are complex valued. To us this is an unphysical contribution which 
we neglect. After all, the discontinuity is directly related to the amplitude to produce 
physical particles, albeit with an (imaginary) chemical potential for color. 

By using the decomposition of the spectral functions, Eqs. (84) and (85), 


tr 


P*a{^,k)p{uj-E,k-p) 


47re(a; — E)6{{P — K^) k){u — E — k + k ■ p) 

+p*_^{u, k){u - E + k - k ■ p) 

(93) 


where p = \p\. 

Since k ^ T, hj using the assumption p ':$> T, we hnd p ^ k. Using this and 

= 0, 

e{u - E)6{{K - Pf) ~ - ^5 (cos 0. (94) 

Zi'pK \ K / 
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Thus, 


Disc i7;(P) ~ + 2 gj ^ y j ^ K(-P) - n,(0)) -^5 (cos 6 - ^ 


X k) i^-1 + ^ j + P-a{<^, k) i^-1 - - 

“2 -na(0))y dk k j 


du 


u 


X [Plai^^k) -1 + ^ +p\{u,k) -1 - ^ 


u 


(96) 


where we introduce an ultraviolet cutoff, p. It is useful to use the sum rules [40], 


r t A 1 

r w A -u i. 

J = ■ 

Using the spectral functions in the time-like region, Eq. (86), we obtain 

1 


(96) 

(97) 


Disc nj^(P) ~ + 2 ^9/ ^ (n„(-£)) - n„(0)) 

f 


X dkk 2 


Z+a{k) ( —1 + 


(^+a{k) 

k 


+ Z_a{k) -1 - 


U_a{k) 


2 ^ K(-^) - ^<^(0)) 


X 


p{uj-a{p) - (^+a{p)) + 2 / dk{uj+a{k) - UJ-a{k)) 


(98) 


where we have used [40] 


iuj±Tk){ul k"^) 2 -^± k 

m^k dk 

(99) 

The wave function constants and the mass shells are functions of the color index, a. 

but we suppress this index for now to make it easier to read. 

form for the mass shells at hard momenta, k 3> gT, 

By using the asymptotic 

n;+ ~ A: + ^ , 

(100) 

uj_ cx k , 

(101) 
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we get 


Disc77;(P) ~ + 2'^gj - n^{0)) 

( 102 ) 


We now make the dependence of the thermal qnark mass on the color index a manifest 
again. 

Evalnating the integral by nsing Eqs. (100) and (101) at the leading-log accn- 

racy, 

Disci7/)(P) ~ + 2 ^ ^ ^ ml^a (na(-P) - n„(0)) In • (103) 

The lower limit of the integral comes from k ~ gT, in which Eqs. (100) and (101) 
becomes nnreliable. 

To leading logarithmic order, then, the Q“’s only enter throngh the statistical 
distribntion fnnctions of the qnarks, and the qnark thermal mass. By nsing Eqs. (24) 
and (45), the contribntion to the prodnction rate for photons from soft qnarks is 
fonnd to be 


a=l 


X 


■1+2 Tdk 


m 


qk^ 


dr rsoit(n\ dr 


pQGP 


The resnlt in the pertnrbative QGP [4] is 

dr 


P 


d^p 


— ^ Qf o aemOs 2 

pQGP f ° ^ 


n{—p) — 1/2 
1 - 


In 




g2T2 


(104) 


In the semi-QGP, this is modihed by a Q-dependent factor, 

1 Ea"^qka(^a(-P) -na( 0 )) 


fT{Q) = 


N ir{—p) — n(0)) 


(105) 


where mqk is the thermal mass when = 0, whose expression will be written later. 

To evalnate the photon production rate in the semi-QGP, we need the explicit 
form of the thermal quark mass when Q°^ ^ 0. From Ref. [21], 

2 / ^ ~ 1 _ \ 

” 24 (^ ■ r - -^(Q*)) - jj (-^(o) - AQn) j ■ ( 106 ) 
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The function A{Q) is given by 


A{Q) 


dEE 


TT^ 


+ 


^(E+iQ)/T — ^{E—iQ)/T — 


( 107 ) 


and A{Q) = A{Q + nT). Note that A{Q) is an even function of Q. 

Our dehnition of A{Q) differs by E from that in Ref. [21], which we do to 
emphasize the physics in the following section, Eqs. (122) and (123). Also for the 
purposes of this discussion to follow, we note that in Eq. (106) the terms involving 
Eb A{Q‘^ — Q^) and ^(0) are from the gluon distribution functions, while A{Q^) and 
A{Q°') are from the quark distribution functions. 

In the perturbative QGP, the thermal quark mass squared is 


m 


qk 


24 


= — A^- — T 


N 




2N 


1 \ 


( 108 ) 


In the hrst expression the 1 is from the gluon distribution functions, while the +1/2 
is from the quark distribution functions. 

It is direct to evaluate A{Q°‘) in terms of the dimensionless variable = 
QV(27rT), Eq. (6), 


+l(Q) = (l-6|gUodi(l-k|modi))T2 . (109) 

While nominally a quadratic polynomial in q, some care must be taken in using this 
expression. Only the absolute value of q enters because by construction Eq. (107) is 
even in Q. Secondly, q is dehned modulo one, since only exp(±27rig“) enters into the 
Bose-Einstein distribution functions in Eq. (107), so the are manifestly periodic 
variables. 

Equation (105) can be simplihed for large photon energy p S> T. In this case, 
na{—p) ~ n{—p) ~ 1, independent of q. We further make use of the fact that all Q’s 
pair up as in Eq. (8) and the corresponding thermal quark masses are identical for 
the components in the pair, = fn^^N+i-a- Consequently, we have 

N N/2 N 

(^a(-p) - na{0)) ~ Y (2 - ^a(O) - n7V+l-a(0)) = Y "^Ska (1 “ ^(0)) . 

a=l a=l a=l 

( 110 ) 
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This allows us to express as the ratio of avergae thermal quark mass sqaured 

when 7 ^ 0 to the perturbative thermal quark mass: 

ff\Q) ^ 


1 Eg^^gka 

N rnl^ 


( 111 ) 


We note that Eq. (Ill) is derived assuming an even N. The conlusion holds for odd 
N also. 

For three colors, taking the eigenvalues as in Eq. (11), the components of ther¬ 
mal mass read 


= ^1^3 = 


= 


6 

9^ 

6 


1 - ^ g + 5 


(l-3g) . 


The suppression factor is then 


fr(Q) = 3 


(1 _ 3,) + I 1 - H , + 5 7 '' 


{n{—p) — 77.(0)) 


For large energy, we obtain a simple polynomial in q, 


fT\Q) ^ 3 


(l-3g) + 2(l-^g + 5g2 


3 3 10 2 

= 1 - 4 (J + y 7 . 


( 112 ) 


(113) 


which agrees with the suppression factor for the hard contribution, /^(Q). Altogether, 
the photon production rate from soft momentum exchange is 


dr 2 1 -p/T 1 

n “ema.4 — e In 


d^p 


f 


TT^ 






f,{Q) • 


(114) 


Comparing the hard contribution in Eq. (79) to the soft contribution in 
Eq. (114), we see that the dependence upon the momentum cutoff p cancels. This is 
a nontrivial check of our computation. The sum of the two contributions is 


dr , . dr 

’’IFp = 


pQGP 


where 


p 


dr 

d^p 


pQGP 


/ 


P 

g^T 


(115) 


(116) 


We can extract Q from lattice results of Polyakov loop and obtain f{Q) as a function 
of the temperature. The result is shown in Fig. (7) . 
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f(Q) 



T/MeV 


FIG. 7. The suppression factor f-y{Q) versus temperature, with the loop from Ref. [26]. 

E. Why so few photons are produced in the semi-QGP 

For dilepton production we found a moderate enhancement near T^. In contrast, 
Fig. (7) shows that photon production is strongly suppressed in the semi-QGP, versus 
the perturbative QGP. To understand the suppression of photons, as in Sec. Ill B it 
helps to generalize the computation to an arbitrary number of colors. In the calcu¬ 
lation of the contribution from the Gompton scattering, the following product of the 
distribution function appears, as was discussed in Sec. IV G: 

1 e-(Ei-iQb)/T^-{E2-iQc+iQt)/T(^^ _ . (117) 

b,c 

Here the Boltzmann approximation was applied to the initial state, and we took the 
large-iV limit, in which we ignore the second term of (T'^'^)be (Eq. (3)) appearing in 
the matrix element squared. The factor 1/iV^ was multiplied for normalization. The 
quantity above becomes 

^ _ n^{E')) (118) 

C 

after partial cancellation of the phase of the distribution functions in the initial state. 
Here we note that this cancellation is not complete unlike the dilepton case: The 
phase iQc/T still remains in the present case while the phase completely cancels for 
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dilepton production in the Boltzmann approximation [29] . By performing the sum as 
in the dilepton case and using Eq. (10), this expression can be rewritten as 

n{NE') (119) 


in the conhned phase. We see that this expression vanishes in the iV —>■ oo limit, 
unlike the dilepton case. The origin of this behavior can be tracked to the fact that 
the cancellation of the phase of the distribution functions for the initial state is only 
partial, and not complete. This is because that the initial state for photon production 
is not a color singlet, as it is for dilepton production. 

For the contribution from pair annihilation, the product of the distribution 
functions is, in the conhned phase and the large-limit, again e~^^'^n{NE'). We 
note that previously [29] , we gave a similar but simpler analysis, using the Boltzmann 
approximation to both the hnal as well as the initial state. 

Next, let us discuss more quantitative point: the origin of the 1/N^ dependence 
of the suppression factor in the conhned phase. For hard photons, with E ^ T, 
we have shown that the ratio of photon production in the semi-QGP, to that in the 
perturbative QGP, is just the ratio of the thermal quark masses squared, of course 
summed over color: 


N 




“i ^ ..1 


m 


qkci 


( 120 ) 


This result is not surprising, as the photon production rate is usually written [4, 5] 
as proportional to the thermal quark mass squared. In the perturbative QGP this 
is somewhat trivial, however, as photon production is naturally proportional to ~ 
e^g^T'^. This relation is less trivial in the semi-QGP, since then the thermal quark 
mass is a function of the Q“’s. Of course Eq. (120) holds only to the order at which 
we compute, which is leading logarithmic order. 

To illustrate how large photon suppression can be, we take the most extreme 
case, the conhned phase. It is most useful to use the form of A{Q) in Eq. (107), as 
an integral over the energy, E. The gluon distribution enters as 

^ . f, 1 „ 1 

( 121 ) 


1 


Aicr - Q^) = r dE E E 

a,b=l “'0 a,6=1 


^{E-i{Q--Qb))/T _ X 
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In the perturbative QGP, = 0, this integral is .4,(0) = T^. In the conhned phase, 
we use Eq. (18) to obtain 

N . R 1 t 2 

( 122 ) 


1 

JP 


6 

- Ocont) = - / dEE 

a,b=l ^ 


^NE/T 


T2 

JP ■ 


Notice that the integral over E is exactly the same as when = 0. The only difference 
is that because only loops which are multiples of N enter, the energy enters not as 
E/T, but as NE/T. Hence in the conhned phase we can replace T by T/N: as the 
integral is ~ T^, this term is suppressed by 1/N^ relative to that in the perturbative 
QGP. 

From Eq. (122), we see that the terms involving the gluon distribution function 
in the thermal quark mass squared, Eq. (106), cancel identically. This leaves only the 
terms from the quark distribution functions, which are functions of the color index a. 
However, photon production only depends only upon the sum over colors, Eq. (120), 
and so we compute 

. N 


ETOL,)"-^/ dEE±Y: 


N 


b=l 


6 

= - — dE E 


N ^ e(®-*QcTnf)/'r + 1 
6=1 

1 T2 


TT^ 


(123) 


eNE/T + 1 2iV2 

by using Eq. (17). Again, this result is precisely 1/A^^ times the result for (5“ = 0. 
We thus hnd that in the conhned phase, the square of the thermal quark masses, 
summed over color, is 

^ 1 \ rr2 


^^mqk(Qconf) — 


N 


a=l 


A^^-l\ 1 

2N ) 24 IP 


(124) 


Gomparing to the thermal quark mass in the perturbative QGP, Eq. (108), we obtain 

1 


MQ 


conf ) — 




(125) 


The coefficient of 1/3 arises as follows. As discussed following Eq. (108), for the 
thermal quark mass squared in the perturbative QGP, the gluon terms contribute 
two thirds (the 1), and the quarks, one third (the 1/2). In the conhned vacuum the 
gluon distributions cancel identically, while the quark terms are precisely 1/A^^ times 
that for Q°‘ = 0, or 1/(3A^^) in all. 
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This shows that photon production is strongly suppressed in the confined phase, 
by 1/A^^. Because the coefficient is small, = 1/3, even for three colors the suppres¬ 
sion is significant, = 1/27. This is why the suppression in Fig. (7) is so dramatic. 


The above analysis applies to the soft contribution to photon production. It 


can also be computed from the hard contribution to photon production, since the 
suppression factor is common. As demonstrated in Sec. IV C, there are two contri¬ 
butions. That from Compton scattering is given in Eq. (74), where by definition, 
/comp(O) = 1 in the perturbative QGP. To compute its value in the confined phase 
of the pure gauge theory, we remember that the only nonzero loops are those which 
wrap around a multiple of N times, Eq. (10). Hence 



The contribution of pair annihilation is given by /pair(Q) in Efi- (78), where /pair(O) = 
1. In the confined phase. 



(127) 


and the contribution from pair annihilation vanishes identically. 

In the confined phase, then, the hard part of photon production only receives 
a contribution from Compton scattering. From Eq. (80), relative to the perturbative 
QGP photon production in the semi-QGP is one third the sum of Compton scattering, 
plus equal contributions from pair annihilation in the t and u channels. Since pair 
annihilation vanishes in the confined phase, Eq. (127), and the contribution from the 
Compton scattering is just 1/N^ times that of the perturbative QGP, Eq. (126), in 
all we obtain a relative suppression factor of 1/{3N^), Eq. (125). 
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FIG. 8. The diagram of four point function with labels indicating fields insertions on 


different branches of the Keldysh contour. 

V. COLLINEAR RATE 

A. Review of AMY’s calculation of photon production 

1. Photon self-energy in RA basis 

First we recapitulate the detailed analysis by Arnold, Moore, and Yaffe (AMY) 
[5] for the collinear contribution to the photon production, in the case that Qa = 0. 
We start with the expression for differential photon emission rate, Eq. (45). In the 
1/2 basis in the real time formalism, is given by 



X Gii22(-i^l, P + K,,-P- A2, A 2 ) . 


( 128 ) 


As will be justihed in the next subsection, (P+2iFi)^ and (P+2A'2)i^ come from quark- 
photon vertices. G'ii 22 (—A'l, P-\-Ki, —P—K 2 , K 2 ) is the Fourier transform of the four- 
point function Gii 22 (a^i, 2 : 2 ,2/i, |/ 2 )- The labels 1, 2 distinguish different held insertions 
on the Keldysh contour. Fig. (8) summarizes the held labeling and momenta how, 
with convention that the momenta how from right to left in propagators. It is easier 
to calculate the four-point function in the RA basis, which is dehned for quarks and 
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gluons as 


'01 + '02 , , , 

Vr= 2 , VA = wi 02 ; 

(129) 

1 4/^ 

/IM _ 1 ^ 2 /iM 

73 r — , ^A — ^2 ■ 

(130) 


In this basis, Gra and Gar correspond to the retarded and advanced propagators, 
respectively. The propagator Gaa vanishes, while vertices with an odd number of 
R indices vanish. The latter is true for quark-gluon vertices only, but only these 
are relevant to the calculation of the collinear rate. To perform the calculation in 
the RA basis, we need to decompose G 1122 in terms of four-point functions. While 
there are in total 16 four-point functions, only 7 of them are independent [41]. The 
decomposition into an independent set has been done for neutral scalar in Ref. [41]. 
It is easily generalized to the case of fermions with /i = 0 as 

Gn22 = chiGaarr + Oi2GAAAR + ocsGaara + CH^GarAA 
+ OibGRAAA + OCqGarrA + OC/GarAR 
+ f^l^AARR + I^^Gaaar + (^3G*aarA + AGarAA 
+ (^^GraAA + f^G^ARRA + (^7^ARAR ■ ( 131 ) 

Detailed analysis by AMY [5] shows that the collinear rate receives contributions only 
from Gaarr- Thus the only coefficients which we need are cti and /9i, 

ai = u(pi)n(p2), A = -(1 - - n{pi)) | ^ ^ . (132) 

-1 + n[p2,) + n[pi) 

In our case {pi,P2,P3,Pi) = {-Ki,P + Ki, -P - K 2 ,K 2 ). 

For the four-point function Gaarr, the RA labeling is uniquely fixed as in 
Fig. (9). The contribution to the collinear regime arises from the kinematic regime 
where hard quark in the loop is nearly collinear with the photon: ~ k\\ > T, 

k± ~ gT^ with || and T defined with respect to photon momentum p. The gluon 
exchanged between the quark lines are soft: gT, q ~ gT. From the collinear 

scattering of quarks and gluons, the energy of the quarks remain unmodified at order 
T. With this kinematic simplihcation, Eq. (132) reduces to 

tti ~/3i ~ u(/c||-Fp)(l - n(fc||)) , (133) 
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FIG. 9. The four point function Gaarr with labels of all internal lines uniquely fixed in 
RA basis. 

and therefore, 


G1122 — 2 tti Re Gaarr ■ 


(134) 


2. Reduction of spinor structure 


To proceed, we need to see how collinear enhancement works. To do so, consider 
the convolntion of two qnark propagators, which enters as a nnit npon inserting an 
additional glnon scattering into Gaarr- 

/ AUO 

—SAR{P + K)SnA{K). ( 135 ) 


Sar and Sra are advanced and retarded dressed qnark propagators: 

1 17° — 7-fc l7° + 7- fc 

— 2 Aq — Ag 2 Aq + Ag 

where 

Ao = , Ag = \k-t\. 


( 136 ) 


( 137 ) 


Here ^ = S^ 7 ^ is the retarded or advanced qnark self-energy, and does not have a 
spinor strnctnre. Note that dne to rotational symmetry, k || S, so Ag = \k — S|. The 
advanced and retarded propagators differ only in the sign of the damping rate, which 
corresponds to the imaginary part of self-energy S. Both Sar and Sra have two 
poles with positive and negative energies. The collinear enhancement occurs when 
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two poles coming from the two propagators pinch the real axis of ko plane. Thus, it 
suffices to consider the pole contribution: 


1 7° —7-fc 


for > 0 


— S(k) ~ ^ 2 Aq-As 2kO{Ao-As) 

' 1 7°+rfc ^ ^ for < n 

It is useful to write ^ in terms of spinor sums: 

jt = '^^Us{k)Us{k) , for > 0 , 

S 

= ''^Vt{k)vt{k) , for < 0 , 

t 

where u and v refer to the spinor basis 

u ^ , 

V ^ (Vk ■ -^k ■ a'qA , 


(138) 


(139) 


(140) 


with 


^^2) : ^12) 1 1 , 2 . 


(141) 


Here, = (1, —a*) with Pauli matrices a*. Note that k^,k '>T and 

gT, so we can take 

Aq - As ^ k^±^Tk - Ek , Ho + ~ + .Efc , Ek = \/kP‘ + m? . (142) 


The asymptotic thermal mass m and damping rate rfc/2 are of order m ~ gT and 
Tfc ~ g^T^ and the explicit expressions of these quantities will be given later. For the 
retarded (advanced) propagator, we take the positive (negative) sign, respectively. It 
is not difficult to hnd that the pinching of poles occurs when kP ~ k \\, with k\\ dehned 
with respect to photon momentum p. With the approximation in Eqs. (142) and 
(135) evaluates to 


dk^ 

—Sar{P + K)Sra{K)c^ 

ZTT 




47 (7 +p) iV + i5E) 


(143) 


k'^=k\\ 


where 

T = ^(Tfc + Tfc+p) ,6E = Ek sgn(/c||) + p - sgn(fc|| + p) . (144) 
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Note that and fC in the nnmerator of Eq. (143) carry independent spinor indices, 

which are to be contracted with qnark-glnon and qnark-photon vertices. Contract¬ 
ing each qnark-glnon vertex with two spinors from the propagators joining it, with 
Eq. (139) there is one of two sitnations, depending on the sign of /c°. Since 

Ut{K)-f^Us{K)=2K^^ 5ts , Vt{K)^^Vs{K) = 2K^^ 5ts, (145) 

each gives the same resnlt. We have neglected the momentnm of the exchanged soft 
glnon Q, since it is negligible compared with K. Cross terms between u and v are 
not allowed becanse mnltiple scatterings with soft glnons do not change the sign of 
k^. We have not inclnded the conpling constant g and color factors, which will be 
discnssed separately in the next snbsection. According to Eq. (145), each qnark-glnon 
vertex gives rises to 2K^, while maintaining the qnark’s chirality. 

Now consider the qnark-photon vertex, contracting the left/right qnark-photon 
vertex with two spinors from the propagators joining them. As an example, consider 
fcO > 0: 


+ P), Us{K + Q + P)Yut{K + Q) . (146) 


Snmming over spinor indices and (transverse) photon polarizations, after some alge¬ 
bra [5] we obtain 


Ut{K)^^Us{K + P)us{K + Q + PWut{K + Q) 

s,t=l,2 i=± 


= 4 (/c° +p) k^- {k±_ g_L) 


/(fc0)2 + (fc0+p)2\ 
V (A;0)2(F+p)2 ) 


~ 4A;_l ■ (fc_L -h g±) 


fc; + (fc||+p)n 
k\\{k\\+p) j 


(147) 


The other cases are similar, with the same result as Eq. (147). Note that by dehnition 
P± = 0, so we can write Eq. (147) as 


E E Ut{K)Yu,{K + P)us{K + Q + P)j%{K + Q) 

,4=1,2 i=± 

= E (2A- + Pn2K + 2Q + py . ( 148 ) 
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Apart from an overall factor (fcjj + {k\\ + p)^)/(fc||(fc|| + p)), Eq. (148) allows us to 
interpret {2K + P)^ and {2K + 2Q + Py as quark-photon verices on the left and right 
of self-energy diagrams. We have thus shown in Eqs. (145) and (148) that quark-gluon 
and quark-photon vertices can be simplihed as 2{Ki + K 2 )^ with Ki and K 2 being 
the incoming and outgoing momenta of quarks. We note that = ~9^lv as 

in Eq. (45). 


3. Color structure in the double line basis 

The color structure of the gluon propagator is given by P^^ as in Eq. (171). 
Thus the color sum which appears when a gluon propagator is sandwiched between 
two quark-gluon vertices is 

(r”),, PS" = (149) 

This can be simplihed further by noting that in the photon self-energy, Eq. (149) is 
sandwiched with the quark-photon vertex in the vertices which are all the way to the 
left or all the way to the right. Starting from the left hand side gives 

^ d/ = ^ - v) d = c-F y . (150) 

where Cp = — 1)/{2N) is the quadratic Casimir for the fundamental represen¬ 

tation. If this is iterated further, each quark-gluon vertex preserves the Kronecker 
delta in color, and generates an additional factor of 6"^. After the last quark-gluon 
scattering, the delta function is color is contracted with the right most quark-photon 
vertex, giving an overall factor of N. 

4- Resummation of infinite self-energy diagrams 

We next resum diagrams with arbitrary quark-gluon scatterings. This is done 
by solving the integral equation illustrated in Fig. (10). The graphical elements are 
the same as those dehned by Arnold, Moore, and Yaffe [5], except that we use the 
double line notation for future convenience. The integral equation shown in Fig. (10) 
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FIG. 10. The diagrammatic equation in terms of graphical elements D, I, M and F. 


becomes 


P) = P) + 


(2ir; 


j M(K, Q, PU,f F(K + Q, O^/K + Q,P). 


(151) 


The color structure can be taken as 


1 


M(K, Q, = — K; M(K, Q, P) , 

F(K + Q, P)et.oi = h F(K + Q,P), 
D^JK,P) = D‘‘(K,P)S^, 


(152) 


SO that Eq. (151) simplifies 
Df^{K,P) = P{K,P) + j 


d!^q 


M{K,Q,P) F{K + Q,P) D^^{K + Q,P) . (153) 


Note the color factor from Eq. (150) cancels the factor of l/(2C'i?) in Eq. (152). 

As discussed previously, we regard the quark-gluon vertices and quark-photon 
vertices as {K + Q)^, with K and Q being the incoming and outgoing momenta. As 
a result, 


I>^{K,P) = {2K + Py , 


(154) 


and the function F{K, P) equals 

F{K, P) = {-zf Gar{P + K) 

-1 \ / 1 


^k\\{k\\ + p) J Vr + i^E 
-1 \ / 1 


4/c||(/c|| + p) J Vr + i6E 


Att 6 (2/c° +p- EfcSgn(/c||) - Ek+pSgn{k\\ + p)) 

27r S{k^ — fc||) . (155) 
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The pinching condition is enforced by the delta function. The rung of the ladder 
equals 


M(K, Q, P) = ig2 Cf (2K + Q + 2^)" {2K + QY G^Y(Q) 

~ 4ig^ Cf k, (fc, + p) P" P" G™((3) . (156) 


with = (l,p). The Ward identity and the fact that K'^ is almost collinear with 
was used to simplify Eq. (156). To further simplify Eq. (153), we dehne 

/ dk^ 

—FiK,P)D>^{K,P), (157) 

which leads us to 

{T + tSE)r{k,p)^{2K+Pr + y^C(q,p)/'‘(k + q,p), (168) 

where 

C(q,p) = s^Cf / ^ 2;r 6(," - ,||)(-iG™(Q)P'‘P‘') . (159) 

The delta function in Eq. (159) again results from + g° —/cy — gy) in F{K + Q,P). 
We can further simplify Eq. (158) using the explicit expression of the damping rate 

r, 

= Q^Cf j 2f S(g° - g||)(-*G««(Q)P'‘P-) . (160) 

As this is independent of k,T = {Tk + Vi.^p)/2 = T^. This allows us to write Eq. (158) 
as 


J SE /'■(k.p) = (2K + PY + j C(q,p) [/"(k + q,p) - r(k,p)] . (161) 

As is clear from Eq. (148), only the transverse components of f^{k,p) are needed, so 
we can project Eq. (161) onto the transverse plane, 

i 6E f±(k,p) = 2k± + J [fx(k + q,p) - f±(k,p)] . (162) 


The last element is to determine the propagator, (Q)P^P^. Since q gT, we 
use the HTL-resummed propagator. 


^ ^ 2T 
-^Gj;«(Q)p>‘p'' = — 
9o 



Q‘‘ - ni;(Q) 


(163) 
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where we have taken Feynman gauge, and used go — 5'||- n^, and IIt are the retarded 
longitudinal and transverse self-energies of the gluon: 


nf(<2) = %M^- 

T 

n?(Q) = ^ 


l-^ln 


2g 


go + g 


^^In 

g2 2g ^ 


% + q 

Qo-Q 


(164) 

(165) 


where the gluon Debye mass is given by 

[j + ^y ( 166 ) 

The Wightman correlator for two electromagnetic currents can be expressed as 
= (-)2iVei ^ Re / 1)1 /-(R, P) F(K, P) P) 

= Ne^Y.‘i‘l j -f^3A(kt,p)’Re[nK,P)nK,P)], ( 167 ) 


where 

_ _ + (^11 +Py\ 

g4(fc||,p) =n(fc||+p)(l-u(fc||)) f 2q^kii+py ) ' 

Note that an overall factor of (^y + (^y +p)^)/(^||(^|| +p)) in Eq- (148) is inserted into 
Eq. (167) along with —1 from the fermion loop. From Eq. (45), is contracted 
with —g^i, to give the collinear rate: 


P 


dr^ _ aemlVE/?/ f d^k 


A{k\\,p) Re[2k_L ■ f±(fc,p)] . 


d^p 47r^ J (27r)^ 

Note there is an additional factor of N for each color of W^‘'. 


(169) 


B. Photon self-energy with nontrivial Polyakov loop 

1. Quark and gluon thermal masses with background color charge 

Now we compute the modihcation of the results in the previous section in the 
presence of a nontrivial Polyakov loop. In this case, quantities like the thermal mass, 
the damping rate, and so on are all dependent on the background color charge. This 
changes the color structure of the self-energy diagrams. 
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The quantities relevant for the problem at hand are the quark asymptotic ther¬ 
mal mass and the resummed gluon propagator. The asymptotic quark thermal mass 
is times the quark thermal mass in Eq. (106): 

ml = 2mjka ■ (170) 


Next we consider the HTL-resummed gluon propagator. The resummed gluon 
propagator consists of bare gluon propagators with arbitrary number of self-energy 
insertions, Fig. (11). The bare gluon propagator in the RA basis of the real time 
formalism, in the presence of background color charge, is proportional to 


^ RA^ab,cd 


^ RR^ab,cd 




'jyab 
cd •> 




(P + ie)2-F' ’ -arm 


v: 


ab 
cd ’ 


(171) 


Here we need to recall that in the analysis in the case of Qa = 0, Bose-Einstein 
enhancement was essentially important for the collinear contribution to be as large 
as the 2 —)■ 2 contribution: for soft gluons with k^ gT, n(P) 1/g. This is no 
longer true in the presence of hard background charge Qa ~ T. The only exception 
is for diagonal gluons, Qa = Qb, where Bose-Einstein enhancement is still operative. 
Therefore, we only need the expression of the diagonal components of the gluon 
propagator. Thus we contract the bare gluon propagator in Eq. (171) with 6ab- As a 
result, the soft diagonal gluon carries only one index. 


SabV:^ 


^ab^cd'Pac 7 


(172) 


where we dehned a color projection operator for one-index gluons. 

Vac = 5ac-^, (173) 

which satishes VabVbc = Vac- In terms of one-index projection operator, the resummed 
gluon propagator is given by the sum of the following terms: 

-^Vab + -^VacV-cd-^Vdb + -^VacV-cd-^V deV-ef -^Vfb + . . . , ( 174 ) 

where the Lorentz indices are suppressed for the time being. Here Tied is given by 
UediQ) = Vcc,dd{Q), with Ilcd,ef{Q) is the gluon self-energy. Each color projection 
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operator is accompanied by a momentum dependent part of the bare propagator 
1/Q^. The color structure of one-index gluon self-energy is 


ffaft ^ab^a ^ ^ab 


(175) 


By restoring the Lorentz indices, the gluon self-energy in the HTL approxima¬ 
tion [22] is given by 


Here 


KIAQ) = . 

0 f dQg 


+ q 


(176) 


(177) 


47r Q.Q±ie 

with Q = (l,q). The sign in the denominator is plus (minus) if the self-energy is 
retarded (advanced). The form of 611^'^ is identical to that in perturbative QGP. The 
Polyakov loop dependence is entirely in the gluon Debye mass: 


■^^-^ab,cd g 


N 




.e=l 


Nf ( 


N 


- ( AiQn - J 1.4(0“) + .4(0“) - ^ ^ .4(0 




N 


e=l 


which for Qa = 0 reduces to Eq. (166), namely 

T _ 1 


M'l,, = ( iiV + 1 P- . 


6 


•yab 


(178) 


(179) 


Equation (178) leads us to 


Fr = j {^A(Qn - lV/.4(0“) j 4n'"'(0) . 


9 


6 


N 


= J I iV.4(0“'’) -N,{A{Q‘)+ .4(0*) 


5^ .4(0') 4n'“'(0) 


(180) 


e=l 


where we have restored the Lorentz indices. 

Formally the two terms in Eq. (175) are of the same order if we regard 6ab ^ 1/N. 
However, we show in Appendix B that a naive large N limit is justihed. This allows 
us to disregard the term proportional to G, so that the gluon Debye mass becomes 


Ma = — 

a 3 


N 


Y,A(Qn - N,A(Q‘ 


e=l 


(181) 
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FIG. 11. The resummed gluon propagator as an infinite series of propagators with arbitrary 
self-energy insertions. Each unfilled circle represents a self-energy insertion 11. The double 
line notation is used. 



FIG. 12. The diagrammatic equation for a color index e, in the presence of nontrivial 
Polyakov loop. 


where Ml is defined as = Ml5ab- When = 2'KTq ^ 0, the explicit form of 

Ml ioi N = Nf = ?) are 



(182) 


2. AMY’s integral equation with one color index 

We next generalize the integral eqnation of AMY for a nontrivial Polyakov 
loop in the limit of a large nnmber of colors. From the discnssion of the previons 
snbsection, all elements of the graph carry one color index, as in Fig. (12). 

The integral eqnation analogons to Eq. (153) is 


Di{K,P) = Ii‘(K,P) + 


I 












Most quantities only need trivial modifications: 


W,P) 

FeiK,P) 


{2K + PY , 

_ ^ Ott 

4k\l{k\\+p) Te + tSEe ^ 


h) ’ 


M,{K,Q,P) = 4iytfc||(fc|| +p)-P'‘F"G™ (0) . 


(184) 


A distinct difference is the color factor Cp m changes to 1/2. This follows from 
enforcing color neutrality on the soft gluon and dropping terms 1/A^ in the gluon self¬ 
energy. Apart from this, the color index e enters T and 6E through quark asymptotic 
thermal mass Eq. (170) and gluon Debye mass Eq. (181). Explicitly, 


5E, 


p 

EfcSgn(A:||) +p- Ek+p sgn(A:|| + p) ~ 2k\\{k\i +p) ^^^ ^ ’ 

— (1 -;/) 1“ _ n«^(Q) - (?=-n«.(Q)) ’ 


with 


nL(0) = 

KeiQ) = -f 


l-^ln(^ 
2g \qo-q 


^^In 

2q 


go + g 
go - g 


Following the case with Q = 0, we dehne 


/ dk^ 

— E,{K,P) Di^iK,P). 


Similarly, 

(r, + i6E,)P{k,p) = (2A- + PY + j ^ Ce(q,p) /"(k + q,p) , 

where 


G.(q,p) = p ■ 

The term proportional to Tg can be written in terms of Cg, 

r. = 2^ 4(«, - ?ii)(-!G™ (0)P'‘P''). 


(186) 

(187) 

(188) 

(189) 

(190) 

(191) 
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FIG. 13. The correspondence between the diagrams for quark damping and for soft gluon 
exchange. The left-hand side is the quark self-energy, whose imaginary part gives the 
damping rate of the quark. The right-hand side is the diagram from gluon exchange, which 
is the cut diagram on the left-hand side. The double line notation is used. 

Physically, this is because quark damping is due to scattering off of soft and diagonal 
gluons. We note that, by using q gT and T/q^, Pe is of order g‘^T. This 

is suppressed by 1/iV compared with P ~ g'^NT, which is the damping rate when 
Q = 0. The diagrams of quark damping are the same as gluon rung M, as illustrated 
in Fig. (13). Note that the Bose-Einstein enhancement Exes the color indices as / = e. 

Therefore we have 

/i‘(k,p) = {2K + pr + j^^ c.(q,p) [/r(k + q,p) - /i‘(k,p)] , (192) 

We again need only an equation for projected fe, 

fe±(k,p) = 2k_L + J C'e(q,p) [fe±(k-k q,p) - fe±(k,p)] . (193) 

The last changes are for Eqs. (131) and (134). Following the derivation of Ref. [41], 
Eq. (131) becomes 

Gn22 = CHiGaARR + Oi2GAAAR + CtS^AARA + CH^GarAA 
+ CtbGRAAA + OCqGarrA + (^iGarAR 
+ / 3 iGaarr + P2GAAAR + (^sGaarA + f^i^ARAA 
+ (^^Graaa + I^gGarra + (^7Garar 1 ( 194 ) 

where A is dehned as complex conjugation together with charge conjugation, i.e. 
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flipping the sign of background color charge. The relevant coefficients are 

tti ~ /3i ~ he(/c|| + p)(l - ne{k\\)) . (195) 

As a result, Eq. (134) becomes 

G 1122 = OtiGAARR{Qa) + OIiG*aaRr{~Q a) ■ (196) 

Note that the background charge enters the integral equation Eq. (193) only through 
Eqs. (170) and (181), which are independent of the sign of Qa, as A is an even 
function. We still have 


1^1122 — 2 ai ^gGaarr^Qo) ■ 


(197) 


Finally, the collinear rate is given by 


d/ly 



(fkl 

(27r)2 


Re[2k_L ■ ie±{k,p)] , 


(198) 


where 

_ _ f ^ + ih +py\ 

Ae(fc||,P) =^e(fc|| +P)(l-ne(fc||)) f ^kj^kw + p)^ J ' 

Note that the factor N in Eq. (169) is replaced by a sum over color index e in Eq. (198). 

To summarize, the collinear rate in the presence of nontrivial Polyakov loop 
is given by Eq. (198), with fej_(fc,p) the solution of Eq. (193). All quantities which 
depend upon the background charge are dehned in Eqs. (170), (181), (185), (186), 
and (190). We note that the Polyakov loop effect enters separately in both the 
longitudinal and the transverse parts, as can be seen from Eqs. (186), (198) and (199). 
In the longitudinal part, the Polyakov loop dependence is reflected in the distribution 
function factor ne{k\\ +p){l — ne{k\\)). In the transverse part, the Polyakov loop effect 
appears in the asymptotic quark thermal mass rUe and the gluon Debye mass Mg. 


C. Photon rate in the collinear regime at large N 

To obtain the photon rate in collinear regime, we need to solve Eq. (193). In 
the limit of large N, the collision term in Eq. (193) is suppressed by 1/N. This can 
be understood as follows. Since 5Ef, ~ ml/k\\ ~ g'^T'^N/k^^, the left-hand side of 
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Eq. (193) is of order g'^T'^N/k\\^e±- On the other hand, by using Eq. (191) the terms 
containing Ce in the right-hand side are of order Egfea ~ Thus, at sufficiently 

large iV, terms in the latter are small compared to the former. 

This allows us to solve Eq. (193) perturbatively. The solutions to zeroth and 
hrst order in terms of are easily obtained. In the argument of all quantities, we 
suppress p but indicate k: 


( 0 ) _ 


2k I 


= 


i5E,{y) ’ 


(1) _ 


d?q 


= 


Ce{q) 


2(k i -|- q i) 


2k I 


( 200 ) 


i(5Ee(k) y (27r)3 [i(5Ee(k-|- q) i5Ee(k) 

Only the solution to hrst order contributes to the photon rate. The relevant combi¬ 
nation is 


d?kj 

W) 


Re|2kx ■ fex(k)] = 4 


(fk I 


d^q 


(27r)2 J (27r)= 


Ce{q) 


k^ 


(kj 


q±j 


[5E,(k)2 


5 E,{k) 5 E,{k + q)_ 

( 201 ) 


We note that, because of this truncation, the LPM effect is suppressed. Using the 
sum rules of Refs. [5, 42], 


dq^dq^^ 


Hq" - g") (-^GJ,^e(Q) =T[^- 


2tt 


ql ql + M^J 


( 202 ) 


Eq. (201) simplihes to 
r d^k I 


(27r) 


Re|2kx ■ fe±(k)] = 2g’ 


2/c|| (/cy-|-p) \ ^ f d?k±_ f d‘^q± 


tm! 


X 


P J J (2vr)2 y (27r)2 gKgi + M^) 
1 / k±- (kx + q±) 


k"] + 777.2 V^i + |k + qP + 


2g‘^T fkiiikii + p) 


(27r)2 

1 


p 


dk^ / dq]_ 


M2 


ql{ql + M^) 


X 


mt 


qj + 2m] 


kl + mlykl + ml 2^^{kl + ql + mlY - {2k^qA_Y) 

(203) 


where in the second line we have performed the two angular integrations, and used 
the formula [5], 



1 

\/ (^i + q\ + 


{2k±q±y 


0 . 


(204) 
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Here k± and q± are of order gT. Nevertheless, since the integrand is convergent in 
both the infrared and the ultraviolet, we can extend the range of the integrations of 
k± and q± to [0, cx)]. The result can be expressed in terms of a dimensionless function 
of the mass ratio Mg/me'- 


(Pk_i 

W) 


Re[2k_L • fe±(k)] = 2g‘^ 


2/c||(/c|| +p) 


P 



(205) 


where 


6^ 




ml 

{k^ + mlY 


In 


mlMl 

{k]_ + mlY 


Ml - 2ml M^jml -kj- - VI) 

2\fA{k‘]_ + ml) Ml{ml — kl_ — M|) + A — {kl_ + ml)y/A 


(206) 


with A = — 2Ml{ml — kjj + {kl_ + ml)‘^. The function Q{Me/me) can be 

determined numerically. 

Consequently, the collinear rate can be expressed as 


P 


dr^ 

d^p 






dki 




/ -^ne{kii+p){l-ne{kii))4g^Tg 



kfi + ih +pf 


pz 


(207) 


The hnal /cy-integral can be done as follows: 
dh (kl + +pf 


ne{k\\ +p) (1 - rieikw)) 


271 \ p 

^+oo+iQ, d/ /2/2 + p72 - 2Ql - Ail Q 


' -OO+iQe 


271 




p 


^ 'n(/+|)n 




(208) 


where I = k\\ + p/2 + iQg. The integrand is exponentially suppressed as Re/ —)■ 
±cx), which allows us to shift the integration contour to the real axis. The following 
integration formulas are useful; 



1 p 

271 eP/'^ — 1 ’ 

1 p(47r^T^ + p^) 
“ ^ 12 (eP/^ - 1) 


(209) 


They can be obtained by integrating / n (/ + |) n (—/ + and Pn (/ + |) n (—/ + 

along the rectangular contour bounded by —oo, oo, oo + 27iTi and —oo + 27iTi. Using 
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Eq. (209), we obtain the collinear rate from Eq. (207) 


dr^ 


Qfp 


P- 




2tt‘T‘ + 2p‘ - 6Qi 

d^p TT^ 3p{eP/'^ — 1) 

For a hard photon, where p ^ the collinear rate is simplihed to 

p-p/r 


dr^ Q^emQ^s f Qf 4:Tp ^^ f iff, 


rrif 


( 210 ) 


( 211 ) 


whose parametric behavior is a Boltzmann factor times a term linear in p. This p- 
dependence is consistent with AMY’s analysis withont the LPM mechanism, and the 
analysis at two-loop order [43]. 

Note that the Polyakov loop only enters throngh the snm ^^^(Me/me). In 
Fig. (14) we show the temperatnre dependence of this fnnction when N = Nf = 3. 
While each individual term Q{M^/me) changes with temperature, especially near T^, 
the sum is remarkably flat, with Ylf,G{Me/me) ~ 3 x 0.015 over a wide range of 
temperature. 

From Eq. (211) the collinear rate is not suppressed in the conhned phase. At 
Erst this is a surprising result, and it is worth discussing in some detail. It happens 
because the soft gluon which is radiated is diagonal in color space, so the quarks in 
the initial and hnal state have the same color indices. The distribution factor which 
appears in Eq. (207) is ne{k\\ + p) (1 — ne(f;||)) = n{p)(ne{k\\) — ne{k\\ + p)). For large 
p > 0, this factor is nonzero only when k\\ -|-p is positive, and k\\ is negative [5, 43]. 
This corresponds to pair annihilation, as illustrated in Fig. (15); the other processes 
correspond to bremsstrahlung, and do not contribute in this limit. 

Since k\\ -|-p is positive, /cy is not only negative, but large. Consequently, as p S> 
T, we can use a Boltzmann approximation for the statistical distribution functions: 

^^ne(fc|| +p) (1 -ne(A;||)) = + P) ^e(-/C||) 

e e 

~ Jl Q-{k\\+P-iQe)/T g-(-fc||+iQe)/T (212) 

N ^ 

e 

= , 


Thus the collinear contribution is not suppressed in the conhned phase because the 
phases cancel between the quark and anti-quark. This is exactly the same cancellation 
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§ (MJm,) 



T/MeV 


FIG. 14. Temperature dependence of function ^(Me/mg). In = 3, the background 
color charge is parametrized as Qe = {—Q,0,+Q). Points of different colors in the figure 
correspond to e = 1 (red), e = 2 (blue) and G averaged over three colors (green). 


as found for dilepton production, and rather unlike the color flow for the contribution 
to photon production from 2 —)■ 2 scattering. 


This completes our derivation of photon rate in the semi-QGP, with a nontrivial 
Polyakov loop at large N. The result is a sum of leading logarithmic term from the 
rate for 2 -f-i- 2, Eq. (115), and the collinear rate in the large N limit, Eq. (210). We 
emphasize that rates for 2 -H- 2 and collinear emission depend upon the Polyakov loop 
in completely different ways. When the Polyakov loop is small, the rate for 2 -f-)- 2 is 
suppressed while that for collinear emission is not. We note that our results are valid 
only for small values of and large N. At moderate values of the coupling constant, 
corrections due to the constant under the logarithm become important. At moderate 
values of N, the LPM effect becomes relevant, and will produce cancellations between 
diagrams with different number of loops, suppressing the photon rate in the collinear 
regime. 
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FIG. 15. The collinear pair annihilation process. The curly line with a blob is the gluon 
HTL-resummed propagator. 

D. Why the LPM effect is suppressed in the semi-QGP 

A nontrivial Polyakov loop is understood as from quantum fluctuations in Aq of 
order T/g. This background gauge field affects quarks and gluons in different ways. 
As shown in Ref. [26], it reduces the density of hard quarks. It also acts as a Higgs 
effect for gluons, giving mass of order T to off-diagonal gluons, while leaving diagonal 
gluons massless. The only gluons which scatter off of quarks in the large N limit are 
diagonal, and so are reduced by 1/N. 

In the perturbative QGP, the LPM effect is relevant because the photon forma¬ 
tion time, tp, is comparable to the mean free path. A, of a quark undergoing multiple 
scattering with gluons in the medium. The formation time is the time scale when a 
collinear photon can be well separated from the quark, which is 

1 T 1 , , 

tfp - f213) 

^ 6E kl+m^ gmr ^ ^ 

On the other hand, the mean free path has the same order of magnitude as the 
damping rate of a quark in the thermal bath, with A rsj 1/P rsj l/{g‘^NT), which is 
comparable to tp. 

A nontrivial Polyakov loop modihes the two scales differently. The thermal mass 
m results from interactions of a quark with hard thermal gluons. The Polyakov loop 
suppresses the quark and the gluon density, and thus also m, by a loop dependent 
factor. The damping rate is due to the scattering off of soft gluons, but as these 
are suppressed for ~ T, only the scattering off of soft, diagonal gluons matters. 
Consequently, A ~ l/{g‘^T) times a loop dependent factor, so at large A^, A S> tp. This 
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implies that quarks rarely scatter more than once during the emission of a photon, 
and thus the LPM effect can be ignored. 


VI. SUMMARY AND CONCLUDING REMARKS 

We calculated the production rates of the dilepton and the real photon in a ma¬ 
trix model of the semi-QGP. The main results of this paper are Eqs. (34), (115), and 
(211). The dilepton production rate was found to be slightly enhanced in the conhned 
phase due to a cancellation in the phases of the statistical distribution functions for 
the quark and anti-quark [29]. By contrast, the photon production rate due to the 
2 —)■ 2 scattering is strongly suppressed for small values of the Polyakov loop, as the 
phases in the distribution functions do not cancel. We showed that the collinear con¬ 
tribution to the photon production is suppressed at large N in the semi-QGP, since 
when the Polyakov loop is small, the Q“’s are large, and off-diagonal gluons do not 
experience Bose-Einstein enhancement. We computed the collinear contribution at 
large N, and found that because of a cancellation of phases, like dilepton production 
it is not suppressed even in the conhned phase. 

These results will modify the theoretical predictions for thermal production 
in heavy ion collisions. Gertainly the production rates for dileptons and photons 
are altered. The elliptic how for these particles are similarly modihed, as the total 
elliptic how is an average over all the phases, from the initial state, to the QGP, to 
hadrons. These ehects were previously discussed in Ref. [29]. However, in that work 
the modihcations of photon production from 2 —)■ 2 scattering and from collinear 
emission were not considered separately. Glearly a more realistic treatment is called 
for. 

In the current analysis, the ehect of the conhnement is taken into account as a 
nontrivial value of the Polyakov loop. It is also interesting to consider the ehect of 
the chiral symmetry restoration as well as conhnement [33, 39]. 

For the future, besides doing a more complete analysis of photon production, 
the most urgent problem is to compute radiative energy loss for light quarks. This is 
closely related to collinear photon emission, and so we expect that near T^, it will be 
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dominated by diagonal gluons for processes in which the color phases cancel. 
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Appendix A: Corrections to Boltzmann approximation to thermal distri¬ 
bution functions 

We will argue that correction to (55) is suppressed by additional exponential. 
We illustrate this in case of Compton scattering. The exact thermal distribution 
factors can be expressed as 



Now the ^/-integral becomes 



(A2) 
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where 


2/0 = 



Ay 


2^ytu{t + M + 4:E{x — E)) 
t + u 


(A3) 


and Iq{z) is the modihed Bessel function of the hrst kind. Note that the leading log¬ 
arithmic contribution comes from t ~ s ^ pT, which implies Ay ~ p. Therefore 
we may set /o((m — n)l3Ay/2) = 1. This leads to the following x-integral 



(A4) 


with 6 = (3{m + n)/2 (3{m — n){t — u)/{2{t -|- u)) > 0 unless m = n = 0. Here, 

F(a, b, c, z) is the hypergeometric function. For non-vanishing 6, there is an additional 


exponential suppression factor e Therefore, we conclude any terms with 


non-vanishing m or n is negligible in Eq. (Al), leaving only the term with m = n = 0, 
which corresponds to the Boltzmann approximation. 

Appendix B: Thermal gluon mass in the presence of Polyakov loop 

We regard Nf ~ A^ as a large number. Naively, the F and G terms give the 
same order contribution in N because 6ab ~ however, as we show below, the 

G term is suppressed by 1/N compared to the F term. By plugging Eq. (175) into 
Eq. (174), we obtain the hrst few terms explicitly: 



(Bl) 
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Here Li{L2) are complicated functions linear(quadratic) in F or G. By induction, we 
can obtain the form of propagator with n self-energy insertions 


(B2) 

It is easy to see the 5ab term can be summed as a geometric series, while the 1/N term 
is not summable in simple manner. In any case, the resummed gluon propagator has 
the following color structure 

5abAa{Q) - , (B3) 


with Aa{Q) and Bab{Q) of the same order in N. Note Aa{Q) = — Fa) ^6ab is 

entirely from F, while Bab{Q) has contribution from both F and G. 

Now we insert the resummed propagator into the graphical element M. Focus¬ 
ing again on the color structure, we obtain after summing over gluon color indices: 


- -B.,] 


— ^ef^hg„ 


^fgAf -j;^iAf + Ag) + — J2Ac- -^B 




'^{Bfc -h Beg) ^ B, 




cd 


cd 


(B4) 


In the above, we have suppressed the P dependence of A and B for notational sim¬ 
plicity. Formally all terms are of the same order if we regard 6fg ~ 1/N and sum 
as ~ N. However, we have learned from the case without background color charge 
that the structure of M is ultimately contracted with dgh on the left, which brings 
Eq. (B4) into the following form 




^fgAf 






1 

Jp 


'^{Bfc + Beg) ^3 ^ Bed 


cd 


(B5) 


We see the 6fg in the first term of the bracket becomes redundant. We can replace 
it by 1. Consequenctly, all other terms are suppressed by 1/N. We will keep only 
the first term, which is fortunately easy to calculate. This approximation amounts to 
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dropping the 1/N term in Eq. (178), leading to the gluon Debye mass, Eq. (181). 
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